arXiv:1509.01515v2 [astro-ph.CO] 12 Feb 2016 


Eur. Phys. J. C manuscript No. 

(will be inserted by the editor) 


Results on light dark matter particles with a low-threshold CRESST-II 
detector 


G. Angloher \ A. Bento C. Bucci L. Canonica X. Defay A. Erb F. v. Feilitzsch 
N. Ferreiro lachellini \ P. Gorla A. Gutlein D. Hauff \ J. Jochum M. Kiefer \ 

H. Kluck H. Kraus J. C. Lanfranchi J. Loebell A. Munster C. Pagliarone 
F. Petricca \ W. Potzel F. Probst \ F. Reindl K. Schaffner J. Schieck 

S. Schonert W. Seidel \ L. Stodolsky \ C. Strandhagen R. Strauss \ A. Tanzke \ 

H.H. Trinh Thi C. Turkoglu ^ M. Uffinger \ A. Ulrich I. Usherov \ S. Wawoczny 
M. Willers \ M. Wustrich \ A. Zoller ^ ^ 

^Max-Planck-Institut fur Physik, D-80805 Munchen, Germany 
^Departamento de Fisica, Universidade de Coimbra, P3004 516 Coimbra, Portugal 
^INFN, Laboratori Nazionali del Gran Sasso, 1-67010 Assergi, Italy 

^Physik-Department and Excellence Cluster Universe, Technische Universitat Munchen, D-85748 Garching, Germany 

^Walther-MeiBner-Institut fur Tieftemperaturforschung, D-85748 Garching, Germany 

^Institut fur Hochenergiephysik der Osterreichischen Akademie der Wissenschaften, A-1050 Wien, Austria 

and Atominstitut, Vienna University of Technology, A-1020 Wien, Austria 

^Eberhard-Karls-Universitat Tubingen, D-72076 Tubingen, Germany 

^Department of Physics, University of Oxford, Oxford 0X1 3RH, United Kingdom 

September 04, 2015 


Abstract The CRESST-II experiment uses cryogenic detec¬ 
tors to search for nuclear recoil events induced by the elastic 
scattering of dark matter particles in CaW 04 crystals. Given 
the low energy threshold of our detectors in combination 
with light target nuclei, low mass dark matter particles can 
be probed with high sensitivity. In this letter we present the 
results from data of a single detector module corresponding 
to 52 kg live days. A blind analysis is carried out. With an 
energy threshold for nuclear recoils of 307 eV we substan¬ 
tially enhance the sensitivity for light dark matter. Thereby, 
we extend the reach of direct dark matter experiments to the 
sub-GeV/c^ region and demonstrate that the energy thresh¬ 
old is the key parameter in the search for low mass dark 
matter particles. 

1 Introduction 

Today, there is overwhelming evidence for the existence of 
dark matter. Precision measurements, e.g. of the cosmic mi¬ 
crowave background, ascribe roughly one fourth of the en¬ 
ergy density of the universe to dark matter, five times more 
than ordinary matter [1]. 

However, physicists did not yet succeed to reveal its un¬ 
derlying nature. One of the most favored solutions is the 
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existence of weakly interacting massive particles (WIMPs) 
thermally produced in the early universe. Those particles are 
beyond the Standard Model of particle physics and would 
have a mass in the range of ^ 10 GeV/c^ to 1 TeV/c^ and an 
interaction cross section of the weak scale. The Lee-Wein¬ 
berg bound excludes WIMPs lighter than ^ 2 GeV/c^, since 
they would lead to an overclosure of the universe [2] . How¬ 
ever, in the last years light dark matter particles below the 
WIMP scale gained more and more interest. In particular, 
asymmetric dark matter models are considered a viable op¬ 
tion as they do not only provide an explanation for dark mat¬ 
ter, but also connect its generation to the baryon asymme¬ 
try [3]. 

The Cryogenic Rare Event Search with Superconduct¬ 
ing Thermometers - GRES ST - is one of numerous direct 
searches around the globe aiming to detect dark matter par¬ 
ticles elastically and coherently scattering off target nuclei. 
To shield against cosmic radiation the experiment is located 
in the ENGS (Laboratori Nazionali del Gran Sasso) under¬ 
ground laboratory in central Italy. In CRESST-II scintillat¬ 
ing CaW 04 single crystals are cooled to mK temperatures 
[4, 5]. While most of the energy deposited in a particle in¬ 
teraction induces a phonon (heat) signal, a small fraction 
(< 5 %) is emitted as scintillation light. The phonon signal 
yields a precise energy measurement, whereas the simulta¬ 
neously recorded light signal allows for particle identifica¬ 
tion, a powerful tool to discriminate electron recoils (dom¬ 
inant background) from the sought-for nuclear recoils. The 






2 


signals of phonon and light detector are measured with inde¬ 
pendent transition edge sensors (TES) read out by SQUIDs 
[4, 5]. We call the ensemble of phonon and corresponding 
light detector a detector module. Each detector is equipped 
with a heater to stabilize the temperature in its operating 
point in the transition between normal and superconduct¬ 
ing state. Via the heater we also inject pulses which are used 
for the energy calibration and for the determination of the 
energy threshold. 

This technology yields accurately measurable sub-keV 
energy thresholds and a high-precision energy reconstruc¬ 
tion. Combining these features with the light target nuclei, 
CRESST-II detectors have a unique potential to explore the 
low-mass regime with dark matter particle masses of below 

1 GeV/c^. 

2 CRESST-II Phase 2 

In August 2015 CRESST-II finished its second extensive 
data taking campaign, referred to as phase 2, featuring 18 
detector modules with a total mass of 5 kg. Most of the data 
acquired in two years of measurement time is dark mat¬ 
ter data accompanied by calibrations with 122keV y-rays 
(^^Co-source), high-energy y-rays (^^^Th-source) and neu¬ 
trons (AmBe-source). In 2014 we reported on first results 
from phase 2, analyzing the detector module with the best 
overall performance in terms of background level, trigger 
threshold and background rejection [6]. This non-blind anal¬ 
ysis proved for the first time that CRESST-II detectors pro¬ 
vide reliable data for recoil energies down to the threshold 
of 0.6 keV [6]. Afterwards, we lowered the trigger thresh¬ 
olds of several detectors, achieving 0.4 keV for the module 
used in [6] and 0.3 keV for the module Lise, which is the 
lowest value obtained in phase 2. 

In contrast to the module used for the 2014 result [6], 
Lise is not equipped with the upgraded crystal holding de¬ 
sign which uses CaW 04 -sticks instead of metal clamps. The 
upgraded holding scheme avoids any non-scintillating sur¬ 
face in line-of-sight to the crystal, thus efficiently vetoing 
backgrounds induced by a-decays on or slightly below sur¬ 
faces [7] . We want to emphasize that the lower threshold of 
Lise is not connected to the holding concept, but solely aris¬ 
ing from a superior performance of the phonon detector. 

Based on the method already discussed in [6] we analyze 
52.2 kg days of data taken with the module Lise (300 g) and 
its threshold set at 0.3 keV. A blind analysis is carried out 
by first defining a statistically insignificant part of the data 
set as a training set. In the previous phase 1 [8] we proved 
that CRESST-II detectors can be operated stably over long 
time periods, in particular concerning the energy reconstruc¬ 
tion. This is also confirmed for Lise as will be discussed in 
section 5. Thus, the training set is well suited to develop all 
methods of data preparation and selection. We then apply 


them blindly - without any change - to the final data set. 
Data from the training set is not part of the final analysis. 

Eor cross checks and validation four independent anal¬ 
yses covering the complete process from raw data to final 
result are performed. Based on the training set results one 
analysis is chosen a-priori. However, the results on training 
set and final data of all different analyses are consistent. 


3 Energy Threshold and Signal Survival Probability 

3.1 Energy Threshold 

In the light of the expected exponential rise of the dark mat¬ 
ter particle spectrum towards lower recoil energies a precise 
knowledge of the energy threshold is mandatory. The finite 
baseline noise causes a deviation of the trigger efficiency 
from an ideal step function. We directly measure the trigger 
efficiency by injecting low-energy heater pulses with a shape 
similar to pulses induced by particle interactions. Eor each 
injected energy Einj the fraction of pulses causing a trigger 
is determined as depicted in figure 1 . 



Fig. 1 Fraction of heater pulses triggering, injected with discrete en¬ 
ergies Einj. The error bars indicate the statistical (binomial) uncer¬ 
tainty at the respective energy. The solid red curve is a fit with the 
sum of a scaled error function and a constant pile-up probability 
Fpiie-up (blue dashed line). The latter is independent of the injected en¬ 
ergy since it is caused by random coincidences of the injected heater 
pulses with any other trigger. The fit yields an energy threshold of 
Eth = (307.3 ± 3.6) eV and a width of a = (82.0 ± 4.2) eV. 


The injected energy Einj is calibrated with 122keV y- 
rays and is verified with low-energy background lines (see 
section 4). The solid curve (red) is a fit with the function 
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where erf denotes the Gaussian error function; /(Einj) de¬ 
scribes the probability for a heater pulse with injected en¬ 
ergy Einj to trigger. The pile-up probability PpPe-up (blue 
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dashed line) accounts for random coincidences of the in¬ 
jected heater pulse with any other trigger and is, therefore, 
independent of the injected energy E-^^y The resulting value 
of Ppiie-up = (5.57 ± 1.17) % is compatible with the expec¬ 
tation based on the trigger rate in this module and the re¬ 
spective coincidence time window. The error function de¬ 
scribes a convolution of a step function for an ideal thresh¬ 
old £’th with the baseline noise modelled by a Gaussian. The 
fit yields an energy threshold of £’th = (307.3 ± 3.6) eV and 
a width of (7 = (82.0 ± 4.2) eV. The latter is in good agree¬ 
ment with the widths of low-energy y-lines (see table 1) con¬ 
firming that the resolution at low energies is dominated by 
baseline noise. 


3.2 Signal Survival Probability 

For all cuts applied to the data, energy dependent signal sur¬ 
vival probabilities are determined by performing the cuts on 
a set of artificial pulses with discrete energies. These are 
created by superimposing signal templates and empty base¬ 
lines, i.e. random triggers periodically sampled throughout 
the whole data taking period. We obtain the signal templates 
by averaging pulses of equal energy. For this analysis y- 
events in the 122 keV ^^Co-calibration peak are used. The 
templates of phonon and light detector are scaled respec¬ 
tively to achieve the desired energy. A dedicated neutron 
calibration confirms that pulse shape differences between 
electron and nuclear recoils are negligible. Therefore, by ac¬ 
counting for the reduced light output of nuclear recoil events 
(quenching) artificial signal events are created. 



Fig. 2 The solid lines represent the signal survival probability after 
successive application of the cuts, as discussed in the text. The sim¬ 
ulated pulses correspond to nuclear recoil events at discrete energies 
starting from the threshold of 0.3 keV (data points). 

The fraction of signals with a certain simulated energy 
passing each cut yields the respective survival probability. 


In figure 2 we show the cumulative probability after each 
selection criterion which are discussed in the following. 

The first cut removes time periods where the trigger rate 
exceeds a pre-defined limit. Such periods are mainly caused 
by an enhanced baseline noise contribution from mechanical 
vibrations induced by the cryogenic facility. For technical 
reasons the rate cut is performed together with the stabil¬ 
ity cut, which only accepts time periods with both detectors 
of a module stably running at their corresponding operating 
points. Rate and stability cut in total remove 2.32 % of the 
collected exposure. 

Due to the low interaction probability multiple scatter¬ 
ings are not expected for dark matter particles. Therefore, all 
events coincident with a signal in any other detector module 
and/or the muon veto are rejected by the coincidence cut, 
reducing the exposure by another 1.12%. 

We reconstruct the energy by fitting the signal template 
to the measured pulse of the corresponding detector and 
quantify the fit quality by the Root Mean Square (RMS) de¬ 
viation. Events where a correct energy reconstruction is not 
guaranteed are mainly removed by a cut on the RMS. The 
line labeled quality cuts in figure 2 shows the combined im¬ 
pact of the RMS cut and additional cuts addressing specific 
artifacts (e.g. pile-up events and SQUID resets). 

For Lise the TES of the phonon detector is evaporated 
onto a small CaW 04 crystal, called TES-carrier, which is in 
turn glued onto the main absorber crystal (composite design, 
details can be found in [9]). Due to the different thermal cou¬ 
plings, pulses in the carrier exhibit a different pulse shape 
compared to those in the main absorber, mainly visible in 
faster rise and decay times. In Lise only a few TES-carrier- 
like events were observed in the training set. Nonetheless, a 
strict cut on the rise-time parameter was applied to remove 
carrier events with a very high efficiency. This has some cost 
on the signal survival probability. One of the cross-check 
analyses uses an artificial neural network to tag the TES- 
carrier like events, exploiting further pulse shape parame¬ 
ters. For Lise no significant impact on the result is observed, 
confirming the performance of the cut on the rise time. In 
all cases a non-negligible fraction of potential signal events, 
especially at low energies (< 2 keV in figure 2), is removed. 
Therefore, for CRESST-III^ we will directly evaporate the 
TES on the absorber crystal, avoiding the TES-carrier. 

For Lise we determine a final signal survival probabil¬ 
ity at threshold of 12 % which is roughly two times higher 
compared to the 2014 result [6]. 


^Due to the major detector upgrades foreseen, the experiment will be 
renamed as CRESST-III. 
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4 Energy Calibration and Discussion 

The energy calibration uses heater pulses of several defined 
energies which are injected periodically throughout the en¬ 
tire data taking period [4, 5]. These pulses are of similar 
shape compared to those induced by particle interactions. 
The heater pulses are used to define a linear response over 
the entire dynamic range of the detector. The absolute en¬ 
ergy scale is finally set by evaluating 122keV y-rays during 
dedicated calibrations with a Co-source. 

The energy spectrum after cuts is shown in figure 3. The 
peak at 2.7 keV is due to cosmogenic activation of tung¬ 
sten, the copper fluorescence line can be seen at S.lkeV. 
The dominant double peak around 6 keV originates from an 
accidental illumination with an ^^Fe X-ray source installed 
to calibrate the light detector of a close-by detector module. 
Despite the significant count rate in the double peak, it only 
marginally impacts the sensitivity for low-mass dark matter 
particles because of its high energy and distinct line shape. 

Depicted in figure 3 is the result from a fit of a data- 
driven background model. The model includes the signal 
survival probability, as discussed previously, a linearly de¬ 
creasing background (dashed black line) and the observed 
low-energy lines. The positions of the peaks are free param¬ 
eters in the unbinned likelihood fit of the model, whereas 
the widths of the peaks are described by a linear function 
allowing for an energy dependence of the resolution. 

As mentioned earlier, the artificial nuclear recoil events, 
as used to determine the signal survival probability, are cre¬ 
ated by superimposing baseline noise samples with signal 
templates. For those artificial events the baseline noise - 
equivalent to the resolution at zero energy - is the only ef¬ 
fect causing a finite resolution of the reconstructed energy. 
This method provides a more precise and - due to the high 
number of empty baselines - statistically superior estimate 
compared to a fit of lines in the low-energy spectrum. The 
determined value of Oq =(62 ± 1) eV is included as a fixed 
parameter in the background model, whereas the linear term 
may be varied by the fit. 

A comparison between literature values But and fitted 
energies Efit of low-energy lines can be found in table 1 . Al¬ 
though we do not use these lines for calibration we find an 
excellent agreement with the nominal energies, demonstrat¬ 
ing the outstanding accuracy of our energy reconstruction. 
Listed in the last column of table 1 is the energy resolution 
evaluated at the respective position of the line^. 

For every particle interaction the total energy deposited 
in the crystal is shared between the phonon and the light 
detector. This split-up differs for different event classes, in 
particular for y-rays (used for calibration) compared to nu¬ 
clear recoils (potential dark matter signal). Nuclear recoils 


Table 1 Comparison between literature values and fitted values for the 
low-energy lines observed in the spectrum of the module Lise (fig¬ 
ure 3). Literature values are taken from [10] (weighted with intensity 
if needed). Energy resolutions evaluated at Efit are listed in addition. 


Origin 

Eiit (keV) 

Eflt (keV) 

Efit Efit 

r^iit 

a(eV) 

i79Ta 

2.601 

2.687 ± 0.020 

-h3.3 

79 ± 1 

Mn Ka 

5.895 

5.972 ± 0.002 

+1.3 

101 ± 1 

Mn K« 

6.490 

6.562 ± 0.003 

+1.1 

105 ±2 

Cu 

8.041 

8.133 ±0.034 

+1.1 

115 ±2 


produce less light and, thus, more energy remains in the 
phonon detector. In [6] we successfully demonstrated that a 
precise measurement of the scintillation light allows to take 
into account this effect accurately. Typical light detectors 
operated in phase 2 (like the one used in [6]) have an en¬ 
ergy resolution which is one order of magnitude better than 
the one of Lise. For the energies of interest, the light signal 
in Lise is, therefore, dominated by statistical fluctuations of 
the baseline noise. Thus, we do not apply the correction as 
used in [6] . This causes the slight overestimation of the re¬ 
constructed energies in the order of a few percent, as can 
be seen in column four of table 1 . In any case, we find this 
deviation to have negligible infiuence on the final result. 



Fig. 3 Spectrum of all events (blue) up to 9keV, truncated at a 
bin content of 100. The red line corresponds to a data-driven back¬ 
ground model including a linearly decreasing background (dashed 
black line) and contributions from ^^^Ta (Ml) (2.7 keV), ^^Ee (6.0 keV 
& 6.6keV), and Cu fiuorescence (S.lkeV). The drop towards lower 
energies results from the declining signal survival probability (see fig¬ 
ure 2). This plot illustrates that the detector Lise exhibits an almost fiat 
background level down to the threshold energy. 


5 Long-Term Stability 


^The values obtained are compatible with uncorrelated fits of the lines, sensitivity for light dark matter critically depends on 

not constraining the Gaussian width. the ability to trigger and reconstruct events of very low en- 
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Live Days 


Fig. 4 Trigger efficiency (black, left y-axis) as a function of the com¬ 
plete measurement time. Each data point is determined as the fraction 
of 0.4 keV heater pulses triggering in the respective time bin. For the 
same time bins also the fitted peak positions and the fit errors for Mn 
Ka (blue) and (magenta) from the ^^Fe-source are drawn (scale 
on the right y-axis). Only small deviations from the fitted total mean 
values (dashed lines) are observed. The constant spread of the corre¬ 
sponding data (grey points) illustrates a stable line width. 

ergies. To prove a stable trigger efficiency over time we pe¬ 
riodically inject low-energy heater pulses (corresponding to 
^ 0.4 keV). In figure 4 the fraction of these pulses triggering 
is shown in black together with statistical (binomial) errors. 
The trigger efficiency at 0.4 keV is slightly higher (^95 %) 
in the final data set compared to the one determined in fig¬ 
ure 1 (^90 %). Thus, the value of 307 eV measured for the 
threshold energy is conservative. 

The high statistics in the ^^Fe double-peak allows us to 
monitor the accuracy of our energy reconstruction and the 
stability of the energy resolution at low energies over time. 



Fig. 5 Data taken with the detector module Fise depicted in the light 
yield - energy plane. The solid lines mark the 90 % upper and lower 
boundaries of the e^/y-band (blue), the band for recoils off oxygen 
(red) and off tungsten (green). The dashed blue line corresponds to the 
lower 5 a boundary of the e“ /y-band, events outside are marked with a 
blue circle. The upper boundary of the acceptance region (yellow area) 
is set to the middle of the oxygen band (dashed dotted red), the lower 
one to the 99.5 % lower boundary of the tungsten band. Events therein 
are additionally highlighted in red. 


mentioned before, in this analysis we do not correct the en¬ 
ergy for events exhibiting a different energy sharing between 
phonon and light detector compared to the 122keV y-rays 
used for calibration. Thus, small statistical fluctuations in 
light output cause a slight tilt of the left shoulder of the 
^^^Pb j8-spectrum [6,11]. 


In each time bin a combined fit of two Gaussians is per¬ 
formed to extract the reconstructed peak positions of the 
(blue) and (magenta) lines (see figure 4). The small vari¬ 
ations of the peak positions show the very stable operation of 
the detector. In addition, a scatter plot of the relevant events 
(grey dots) further illustrates that the means and also the line 
widths of the two peaks remain stable over time. 

6 Light Yield 


Nuclear recoils produce less light than electron recoils, 
as quantified by the quenching factor for the respective tar¬ 
get nucleus. The quenching factors, including their energy 
dependence, are precisely known from dedicated external 
measurements [12] and allow to analytically calculate the 
nuclear recoil bands for scatterings off tungsten (solid green), 
calcium (not shown) and oxygen (solid red), starting from 
the e“/y-band. The validity of this approach is, for example, 
shown in [8] and was confirmed with a neutron calibration. 


In figure 5 we present the light yield - defined for every 
event as the ratio of light to phonon signal - as a function 
of energy. Recoils on electrons have the highest light out¬ 
put, set to one by calibration (at 122keV). The solid blue 
lines mark the 90 % upper and lower boundaries of the e“ / y- 
band, with 80 % of electron recoil events expected in-between, 
The band is determined by an unbinned likelihood fit to the 
data [7] . Two prominent features can be seen in the e“ /y- 
band, the double peak at ^ 6 keV originating from the ^^Fe- 
source and a j3-decay spectrum from an intrinsic contam¬ 
ination in the crystal with ^^^Pb starting at 46.5 keV. As 


There are four events (blue circles) in the data well be¬ 
low the 5 (7 lower boundary of the e“/y-band (dashed blue). 
These events are statistically incompatible with leakage and, 
thus, indicate the presence of an additional source of events 
in this module, apart from e“ /y-background. As discussed 
in the introduction, Lise is of a conventional module design 
using non-scintillating clamps to hold the crystal and, there¬ 
fore, is not capable of vetoing all events originating from a- 
decays on or slightly below surfaces [7, 8]. Additional data 
from Lise (before lowering the threshold) and the ten other 
conventional modules further underpins this explanation. 
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7 Leakage into the Acceptance Region 

7.1 Acceptance Region 

Before unblinding we define the acceptance region (yellow 
region in figure 5) to extend from the threshold energy of 
307 eV to 40 keV. In light yield the acceptance region spans 
from the 99.5 % lower boundary of the tungsten band to the 
center of the oxygen band (dashed-dotted red line in figure 
5). For the calculation of the exclusion limit all events inside 
the acceptance region (highlighted in red) are conservatively 
considered as potential signal events. In [6] we found that 
the center of the oxygen band is the optimal choice in terms 
of limiting the leakage from the e“ /y-band on the one hand 
and sampling the signal of all three nuclear recoil bands on 
the other hand. This finding was confirmed for Lise based on 
training set data. However, the impact of the choice of rea¬ 
sonable light yield boundaries is on the level of the expected 
statistical fluctuations. 


7.2 Background Leakage 

In the energy range from 1 to 40 keV we determine an aver¬ 
age background level of 13 counts/(keV kg day). If the con¬ 
tribution of the ^^Fe-source is excluded we obtain a value of 
8.5 counts/(keV kg day) which is still well above the value 
of 3.51 counts/(keV kg day) [13] determined for the module 
used in [6] . These differences almost completely arise from 
different levels of crystal-intrinsic radioactive contamina¬ 
tions. The fact that the crystal of the module Lise is from an 
external supplier while the crystal analyzed in 2014 [6] was 
grown at Technische Universitat Munchen further stresses 
the importance of in-house crystal growth [14]. 



Energy (keV) 


Fig. 6 Energy spectrum of all events in the acceptance region (307 eV- 
40keV, see figure 5) truncated at a bin content of 30 for reasons of 
clarity. For the final result all events are conservatively considered as 
potential signal events to extract an exclusion limit using Yellin’s opti¬ 
mum interval method. 


Since the widths of the bands are dominated by the en¬ 
ergy resolution of the light detector, the modest light detec¬ 
tor performance of the module Lise results in a low discrimi¬ 
nation power between the sought-for nuclear recoils and the 
dominant e“//-background. Thus, we observe leakage of 
e“//-events into the acceptance region, leading to the en¬ 
ergy spectrum of all accepted events depicted in figure 6. 


8 Expected Signal Composition on the Target Nuclei 

To calculate the expected recoil spectrum of dark matter par¬ 
ticles we make use of the following standard assumptions 
on the dark matter halo: Maxwellian velocity distribution, 
asymptotic velocity of 220km/s, galactic escape velocity of 
544km/s and local dark matter density of 0.3 GeV/cm^. We 
use the Helm parametrization of the form factors [15] quan¬ 
tifying effects of the nuclear substructure on the scattering 
cross-section. Finally, the choice of the acceptance region in 
the light yield - energy plane and the energy resolution at 
zero energy (baseline noise) are taken into account. For the 
latter we use the value of op =62 eV as determined by the 
evaluation of the artificial signal events (see section 4). 



Fig. 7 Expected number of counts in the acceptance region for a dark 
matter-nucleon cross-section of Ipb and standard astrophysical as¬ 
sumptions as a function of the mass of the dark matter particle. The 
solid black line depicts the overall count number which is the sum of 
the number of scatterings off the three target nuclei (colored lines, see 
legend). 

Figure 7 shows the number of expected events to be ob¬ 
served by the detector Lise in the final data set as a function 
of the mass of the dark matter particle. The lines are normal¬ 
ized to a dark matter particle-nucleon cross section of 1 pb 
(= 10“^^ cm^). The colored curves (see legend) correspond 
to the three target nuclei oxygen, calcium and tungsten. The 
total sum of the three constituents is drawn in solid black. 
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Due to the anticipated -dependence of the cross sec¬ 
tion, dark matter particles are supposed to dominantly scat¬ 
ter off the heavy tungsten. The energy transferred in the scat¬ 
tering process is a function of the reduced mass of target 
nucleus and dark matter particle. Thus, for a given mass of 
the dark matter particle the fraction of the expected energy 
spectrum above threshold depends on the mass of the target 
nucleus. 

As a result, for dark matter particles with masses above 
5 GeV/c^ recoils off tungsten are expected to be far more nu¬ 
merous compared to oxygen and calcium. For lighter masses 
a substantial part of the tungsten recoils have energies be¬ 
low threshold leading to a strong decrease of the number of 
counts. This results in a mass range completely dominated 
by scatterings off oxygen, because the drop for oxygen and 
calcium is shifted towards lower masses (see figure 7). 

In the limit of very low masses, the reduced mass con¬ 
verges to the mass of the dark matter particles, causing less 
pronounced differences in the shape of the recoil spectra on 
the different target nuclei. This effect is further augmented 
by the influence of the baseline noise. Since the A^-scaling 
of the cross sections still persists, scatterings off tungsten 
account for a slightly larger proportion of the total expected 
signal again. 


9 Result, Discussion and Outlook 

For each dark matter particle mass we use the Yellin op¬ 
timum interval method [16,17] to calculate an upper limit 
with 90 % confidence level on the elastic spin-independent 
interaction cross-section of dark matter particles with nucle¬ 
ons. While this one-dimensional method does not rely on 
any assumption on the background, it exploits differences 
between the measured (see figure 6) and the expected en¬ 
ergy spectrum (see section 8). 

The resulting exclusion limit of this blind analysis is 
drawn in solid red in figure 8. For higher masses this module 
does not have a competitive sensitivity, due to the large num¬ 
ber of background events. In particular, the leakage from the 
^^Fe-source (see figure 6) results in an almost fiat limit for 
masses of 5-30GeV/c^. However, for dark matter particles 
lighter than 1.7 GeV/c^ we explore new regions of parameter 
space. 

The improvement compared to the 2014 result [6] (red 
dashed line) is a consequence of the almost constant back¬ 
ground level down to the threshold which was reduced from 
603 eV to 307 eV. The lower the mass of the dark matter par¬ 
ticle the more relevant these improvements become. With 
this analysis we explore masses down to 0.5 GeV/c^, a nov¬ 
elty in the field of direct dark matter searches. 

The transition point of the dominant scattering target nu¬ 
cleus manifests itself as kink in the corresponding exclusion 



Dark Matter Particle Mass (GeV/c^) 


Fig. 8 Parameter space for elastic spin-independent dark matter- 
nucleon scattering. The result from this blind analysis is drawn in solid 
red together with the expected sensitivity (la confidence level (C.L.)) 
from the data-driven background-only model (light red band). The re¬ 
maining red lines correspond to previous CRESST-II limits [6,18]. The 
favored parameter space reported by CRESST-II phase 1 [8], CDMS- 
Si [19] and CoGeNT [20] are drawn as shaded regions. Eor com¬ 
parison, exclusion limits (90 % C.L.) of the liquid noble gas experi¬ 
ments [21-23] are depicted in blue, from germanium and silicon based 
experiments in green and black [24-28] . In the gray area coherent neu¬ 
trino nucleus scattering, dominantly from solar neutrinos, will be an 
irreducible background for a CaW 04 -based dark matter search experi¬ 
ment [29]. 

curve. Due to the lower threshold Lise starts to be domi¬ 
nated by scatterings off tungsten already at 3 GeV/c^ (see 
figure 7) compared to ~4.5 GeV/c^ for the 2014 result [6]. 

Due to the rather large number of leakage events into the 
acceptance region the result is already not limited by expo¬ 
sure any more. Consequently, only small statistical fiuctua- 
tions are expected. This is confirmed by calculating limits 
for 10,000 Monte Carlo sets sampled from the data-driven 
background model discussed in section 4. The resulting 1 a 
contour is shaded in light red in figure 8. 

In CRESST-III we will substantially size down the ab¬ 
sorber crystals in order to achieve lower energy thresholds. 
Furthermore, we expect two beneficial effects on the light 
signals: Firstly more light reaches the light detector and sec¬ 
ondly the light detector can also be scaled down which leads 
to an enhanced energy resolution. Both improvements will 
increase the background discrimination power. All modules 
will feature an upgraded holding scheme and will mainly 
be equipped with absorber crystals produced in-house due 
to their significantly lower level of intrinsic radioactive con¬ 
taminations. Combining these measures with the enhanced 
discrimination power, a drastically reduced background leak¬ 
age is expected. 

In this letter we prove that a low energy threshold is 
the key requirement to achieve sensitivity to dark matter 
particles of ^(1 GeV/c^) and below. We expect significant 
progress exploring the low mass regime with the upcoming 















CRESST-III experiment, featuring next-generation detectors 
optimized towards the detection of recoil energies as small 
as lOOeV. 
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